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Upper elevations of the higher mountains of the southwestern 
United States are inhabited by several species of boreal mammals 
that have been isolated on these mountain refugia for varying 
periods of time. These restricted populations provide an oppor¬ 
tunity for investigating zoogeographic and evolutionary princi¬ 
ples. One such mammal, the Mexican vole (Microtus mexicanus), 
chosen as the basis of this study, presently inhabits the Transition 
Zone on scattered mountain ranges from southwestern Colorado 
to the Mexican state of Veracruz. It typically is found in diy 
bunchgrass meadows and mesic grasslands scattered among yellow 
pint* and lit at elevations generally above 2000 meters (rn.). M. 
jnexicanus is mote tolerant of xeric conditions than are other mi- 
crotines inhabiting this area and is often found < onsiderable dis¬ 
tances from any source of permanent water (Findley and Jones, 
1962). 1'lie isolated nature of mountains in the Southwest and rel¬ 
atively small geographic ranges of the species have lesulted in a 
series of small, isolated Microtus rnexuanus populations. 

This study concerns four populations of the Mexican vole 
located in the San Mateo, Manzano, Sacramento, and Guadalupe 
mountain ranges of New Mexico and Texas. These mountain 
ranges are separated from one anothei by variable expanses of dry 
lowlands, generally consisting of unsuitable habitat for microtine 
rodents. Three of the four populations have been referred pre¬ 
viously to the subspecies M. m. guadalupensis (Manzano, Sacra- 
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memo, and Guadalupe mountain populations); the remaining 
San Mateo Mountain population was referred to AT rn, mogollo - 
nertsis (Hall and Kelson, 1959). 

I he primaly objectives of this investigation were to reevaluate 
present subs pec Hu assignments of populations of \1. mexicanus 
in the region of study, to compare conclusions obtained by means 
of classical taxonomic methods with those obtained by using 
more recent systematic methods, to examine the /oogeographir 
history and evolutionary relationships of the populations con¬ 
cerned, and to attempt to place the evolutionary c hanges observed 
into a reasonable temporal framework Additionally, this investi¬ 
gation was intended to explore aspects of the biology of M. rnexi- 
canus that previously have not been covered in the literature. 
These four populations of the Mexican vole were examined using 
protein electrophoresis, karyotypic analysis, sperm morphology, 
bacular morphology, and classical morphometries. 

Methods and Materials 

Mexican voles were captured in Sherman live traps during the 
months indicated at the following localities; Texas: Culberson 
Co.: Upper Dog Canyon, Guadalupe Mountains National Park 
(April. August, December); 'The Bowl, Guadalupe Mountains 
National Park (August). New Mexico: Torrance Co.: 4th of July 
Campground, Cibola National Forest (September October); Lin¬ 
coln Co.: South Fork Campground, Lincoln National Forest 
(May, October, November); Soccoro Co.: Bear trap Canyon, Cibola 
National Forest (October). The collecting localities and the extent 
of the mountain ranges involved aie illustrated in Fig. 1. 

After collection, voles were transported from the field to the 
laboratory and were usually sacrificed within two weeks of cap¬ 
ture 1 . A total of 379 spec miens of M. mexicanus were examined, 
including 121 that were live trapped and maintained in the labor¬ 
atory for varying periods before being killed. All were prepared as 
standard museum skins and skulls and deposited in I he Museum. 
Texas lech University. 

Karyotypes were prepared following the in vivo technique of 
Baker (1970). A minimum of five metaphase spreads were counted 
tor each of 42 males and 19 females examined. Terminology des¬ 
cribing chromosomal morphology and fundamental number fol¬ 
lows that of Patton (1907). 

Spermatozoa were* obtained by removing the epididymis from 
freshly killed specimens, mincing it with sc issors, and suspending 
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ranges containing mlloting localities foj populations of .Mirrotus jnexifanus 
examined. Localities are identified h\ imnihcis: 1) Keuitrap Canyon (San Mateo 
Mountains, Souom County, New Mexico); 2) 1th oi July Campground {Mau/auo 
Mountains, Jon.iiui Counts. New Mexico); South Folk Campground (Sacra- 
memo Mountains, Lincoln Counts, New Mexico); 1) Guadalupe Mountains (Cul- 
heison Counts, Texas). Circles indicate approximate collecting legalities within 
each mountain range. 


a small amount of fluid from the minced tissue in an isotonic 
solution of sodium citrate. Several drops of this solution were 
placed on a microscope slide and allowed to ail dry. Slides were 
then fixed in a solution of one part acetic ac id and four parts 
methanol prior to staining with a 0.15 per rent solution of 
Giemsa in hot water. Photomicrographs of spermatozoa were 
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mm 

Eic;. 2. — Batulinn of M/r rotus mr\n anus guadalupensis (IT U 27*326) illwsti :itintv 
tin 1 dimensions described in the text: A-B, Ixuukir length; J)-E, width oi base; 
width oi shall; B-C length oi median distal piocrss; II I, width ot median distal 
pmt ess. 

taken using a Leit/ Wetzlar microscope at a magnification of 
95UX. Measurements in millimeters were taken directly from the 
four by five-inch negatives using Helios dial calipers and a Saku- 
rai map measuring device. Measurements and terminology follow 
those of Linzey and Layne (1974). Ten spermatozoa were mea¬ 
sured for each of the 28 specimens examined, and average values 
were used in the subsequent analyses. 

Penes were removed from freshly killed specimens and stored in 
AFA; bacular preparation followed techniques outlined by Ander¬ 
son (1960). After staining and destaining, penes were dissected 
away from the bacula, and the bacula were drawn in both dorsal 
and lateral views using a camera lucida attached to a Wild M5 
microscope at a magnific ation of 12X. Measurements (in millime¬ 
ters) of 64 specimens were taken directly from the drawings and 
included; bacular length, width at base, width of shaft, length of 
median distal process, and width of median distal process when 
an ossified distal process was present (Fig. 2). 
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1 he morphometric portion of the study was based on .179 sprc- 
imens; only those specimens with a total skull length of 23 mil* 
liineters or more were included in the analysis. This corresponded 
to a minimum total body length of approximately 100 millime¬ 
ters. Four standard external measurements were recorded from 
specimen labe ls and 11 cranial measurements were taken. Cranial 
measurements follow* Coin (1913), Anderson (1951), and Snvder 
(195-1). 

Starch gel electrophoresis was used to assess allozymic variation 
lor 80 Microtus mexicatms fiom the following localities: Upper 
Dog Canyon (11 males, 9 females); South Fork Campground (12 
males, 8 females); 4th of July Campground (8 males, 12 females); 
and Beartrap Canyon (11 males, 9 females). Samples were pre- 
panel from liver, heart, and kidney extracts. Techniques of tissue 
preparation, electrophoresis, and biochemical staining were sim¬ 
ilar to those of Selander et ai (1971), hut staining procedures for 
sorbitol dehydrogenase (SDH) and glutamate dehydrogenase 
(GDI1) were modified after Shaw and Prasad (1970). A summary 
of electrophoretic methods used is given in Table 1. 

Designation of alleles follows that of Smith et al. (1973). The 
allele occurring in the highest frequency at each locus was 
assigned the value of 100 for anodallv migrating systems, or —100 
for those migrating catliodally. Remaining alleles at a locus were 
described as percentages of the 100 allele by comparing relative 
migration distances. When more than one locus was present in a 
system, the most anodal loc us in the system was designated T’ 
and more successively cathodal loci were assigned progressively 
higher numbers. Allo/.yme similarity was assumed if side-by-side 
comparisons failed to establish differences (see Smith et ai , 1973). 

Statistical and clustering analyses were* carried out using SAS-76 
(Ban et at, 1976), TNI VAR (Power, 1970), and NT-SYS pro¬ 
grams. In all analyses, specimens we re grouped into four aggre¬ 
gate populations representing the primary collecting localities; 
these aggregates were considered as Operational Taxonomic Units 
(OTU's). 

SAS programs were used irr univariate analyses to provide 
standard descriptive statistic s and to perform both single classifi¬ 
cation and two-way analyses ol variance in order to test for signif¬ 
icant differences between or among means. When means were 
found to be significantly different, the Sums of Squares Simul¬ 
taneous Test Procedure (SS-STP) w t hs conducted using a UNI VAR 
program to determine 1 maximally nonsignificant subsets. Multi- 
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variate analyses entailed calculation of Pearson product-moment 
correlation coefficients, again using SAS, to compare elertropho- 
letic distance nul correlation matrices with the corresponding 
motphometi ic and sperm morphology matrices, (.luster analyses 
were conducted on the correlation and distance matrices and a 
two-dimensional phetiogram was generated for each using the 
N f SVS program. This program utilizes the unweighted pair- 
group method using arithmetic averages (UPGMA). Phenograms 
were compared with their respective matrices and a coefficient of 
cophenetic correlation was computed for each comparison to 
assess the reliability of the phenogram. A matrix of correlation 
among the characters was computed, and the first principal com¬ 
ponents were extracted. A three-dimensional plot of the OTIJ’s 
onto the first thiee principal components was then prepared. I lie 
percentage of the total variation accounted for by each principal 
component was calculated, as was the contribution of each char¬ 
acter to each of the principal components. Additional data con¬ 
cerning methods and materials may he found in Wilhelm (1977). 

Spec a mkns Kxaminkd 

In addition to specimens collected for this study, other speci¬ 
mens wete borrowed from the* following institutions for examina¬ 
tion; Texas Cooperative Wildlife Collection, Texas A&M Univer¬ 
sity (TCWG); the Museum, Texas lech University (HI;); 
University of New Mexico (UNM); and University of lexas at K! 
Paso (UTEP). Localities are not plotted separately on Fig. 1, but 
are grouped with the four primary collecting sites for purposes of 
analysis and discussion. The list below includes all spec imens 
examined. 

MifTotus mexiranu ;i guatlaiupensts (300).—Texas: Culberson County*. Guada¬ 
lupe Mountains National Park, The Bmx l, 2*1 (6 T I t . 18 TCAVC); Ppper Pojr 
Canyon, 33 (1 J I ); Guadalupe Peak, I (111'), Saw Mkxh o: Otero County: Fini- 
heion, 2 ( l i t 1 ); 8,5 mi. K, 1,5 mi. X Alino Peak. 2 (l NM); 1 mi, N Cloiuktoft, 2 
(FNM); 7 mi. F (Jotuk tolt, 2 (UNM); 7 mi. F, 2.5 mi. N Olomkroh, 7 (l NM); L 
mi S Cloutfunft. 1 (l'i\M): 10 mi. S Clnudrioft. 15 (l NM); 2 mi. \V CJnudrrott, 1 
(lAM); 2.1 mi. \V Cloudcroft, 1 (l I I P); near CInudnoft, 19 (8 l II P. 11 IM P); 
Russian Canyon, 5 mi, S. 2.5 mi. I Chunk iolt, 2 ft I I P): north ol Ruidoso, 3 
(P I FPi. Lincoln County; Padilla Puim, 3 (PTFP): South Fork Campground, 27 
( 1 I P); 5 mi \\ 9 mi. Y Capitate 2 (ONM); 15 mi. W Capiraip 3 (UNM); Capitati 
Mountains, 20 (UNM); Monjtau Peak, 10,000 ft.. 6 (UNM); Lincoln County, ft 
(I I P). Torrance County: 4th ol July Campground, 36 (I IP); 0.5 mi. S Capdlo 
Peak, 9000 ft., 1 (UNM); 5,5 mi. W Tajique, 8 (UNM); Red Gut yon, 0.5 mi. S, 5 
mi. W. Man/iinn, 31 (t NM) Red Canyon, 4 mi. \\\ 1 mi. S Man/ano, 11 (PXM); 
5 mi. W Man/ano. 3 (UNM); l oname County, 4 (TI P). Jf Bernalillo County. Free 
Springs. 8600 li., 21 (UNM). 



['able ].— Electrophoretic techniques utilized in this Uud\. All geh were run for five hours. 
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Microtus mexicanus mogollonensis (79)— New Mexico: Socorro County : Beartrap 
Canyon, 78 (16 TEU, 28 TNM, 4 UI'EP); 1.6 mi. fiom Water Canyon, 1 (UTEP). 


Results and Discussion 
Karyology 

The karyotype of Microtus mexicanus was described from a 
single female specimen by Matthey (1957), who established the 
diploid number as 44. Being unable to identify the sex chromo¬ 
somes, he was uncertain of the fundamental number and listed it 
as a minimum of 54 and perhaps 56. 

The karyotype of Microtus mexicanus given in Fig. 5 shows a 
diploid number of 44 and a fundamental number of 54. No 
chromosomal polymorphism was noted among specimens exam¬ 
ined. The autosomes consist of three pairs of large submetacen- 
trics, one pair of medium submetacentric s, two pairs of small 
metacentrics, one pair of large acroc entries, and 14 pairs of small 
acrocentric s. The X chromosome is a medium-sized submetacen- 
tric, and the \ is a small acrocentric. 

Although no chromosomal polymorphism was found in the 
karyotype of the Mexican vole in this study, two separate chromo¬ 
somal polymorphisms were described in Mexican voles from the 
states of Jalisco and Durango in Mexico (Lee and Elder, 1977), 
and involved differences in both the diploid and fundamental 
numbers. Thus, M. mexicanus is the only exclusively Nearctic 
species of microtine from which infraspecific differences in 
chromosome morphology or number have been reported, 
b Infraspecific variation in chromosome numbei is more com¬ 
mon in Palearctir representatives of the genus, however. Microtus 
hyperboreus and M. middendorfi (Gileva, 1972) and M. juldaschi 
(Bol’shakov et al, 1975) have been reported as polymorphic. 
Other genera of microtines that exhibit chromosomal polymor¬ 
phism include Dicrostonyx torquatus (Rausch and Rausch, 1972; 
Kozlovskii, 1974), Clethrionomys rutilus (Rausch and Rausch, 
1975), and Pitmiys subterraneus (Meylan, 1972). The relatively few 
number of karyotypic studies of Nearctic microtines and the 
apparent uniformity of infraspecific chromosome complements 
might reflect a lack of data rather than extreme chromosomal 
conservatism of the group. 

Sperm Morphology 

Spermatozoa of M. mexicanus from the four populations stud¬ 
ied were similar in all respects. The sperm head, widest just above 
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Fig. 3.— KaryotvfH* of a untie fl I T 273”>t») Mtcrotus tnpxicanits guadalupmsis. 


the base, is asymmetrical with one margin convex and the other 
nearly straight. I he base is smoothly convex with a notch on one 
side. Almost one-half of the head is enveloped by the atrosome, 
which is elongated into a recurved honk that lies on the side of 
the head having the straighter margin. Attachment of the mid- 
piece is to the basal notc h and is therefore eccentric. J he tailpiece 1 
tapers gradually toward the tip and is sometimes difficult to dis¬ 
tinguish from the midpiece. However, with the 1 staining procedure 
used in this study, the tailpiece is generally of uniform appear¬ 
ance, whereas the midpiecc is nearly always granular or mottled 
in appeal ant e. 

Mean values, followed by range in parentheses and sample si/e, 
for selected sperm measurements of Micratus mexicanus follow 
{measurements ne in microns and localities are given in the order 
Bearnap, 1th of July, South Fork, and Guadalupe): length of 
head, 7.77(7.11-8.00) 5, 8.02(7.51-8.51) 5, 8.34(7.60-1). 1 4) 10, 

8.02(8.01-8,78) 8; width of head , 4.86(4.82-1.90) 5, 4.82(4.47-5.05) 5„ 
5.08(1.51-5.05) 10, 5.01(4.81-5.05) 8; length of midpiece , 
18.45(17,12-20.02) 5, 18.10(16.80-18.80) 5. 18.02(16,72-19.06) 10, 
18.46(15.81-20.72) 8; length of tailpiece, 66,02(60.16-71.06) 5, 
70.21(65.81-74.02) 5, 72.02(68.48-76.08) 10. 70.51(65.92-70.01) 8. 

Ihuvariate analysis of sperm data showed a significant differ¬ 
ence 1 in the head width of spermatozoa (/■><(). 05) between some 
localities and in the length of the tailpiece (F<0.01) between oth- 
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cts. The following nonsignificant subsets were generated by SS- 
STP tests: width of head —Guadalupe, 4th of July, Beartrap; 
length of tailpiece —South Fork, Guadalupe, 4th of July; Guada¬ 
lupe, 4th of July, Beartiap. In M. m . mogollonensis , the sperm 
head is slightly shorter but just as broad as that of M. w. guada- 
lupensis. The tailpiece of M. m. mogolloncnsis is considerably 
shor ter than that of the other subspecies. 

The four sperm measurements also were analyzed using the 
NT-SYS multivariate program. Both correlation and distance 
matrices were computed and phonograms representing the phe¬ 
netic relationships were plotted. A Pearson product-moment 
correlation matrix was computed comparing these two matrices 
with the corresponding electrophoretic matrices. Correlation coef¬ 
ficients for neither the distance matrices nor the correlation matri¬ 
ces were significant. The distance phenogram showed that the 
Beartrap and 4th of July specimens had smaller sperm than did 
those from the other two localities. Although on the basis of 
sperm data, the Beartiap (M. m. mogolloncnsis) and the 4th of 
July (AL m. guadalupensis) populations more closely resembled 
each other than either resembled the other two populations, their 
relationship to each other was not of the same magnitude as was 
that of populations from the South Fork and Guadalupe locali¬ 
ties. 

In the principal components analysis, the amount of phenetic 
variation expressed in the first principal component was 65.58; in 
the second, 53.60; and in the third, 0.80. The pel rentage contribu¬ 
tion of each sperm character to each principal component, given 
in the order Component I, II, and III were length of head , 10.51, 
0.71, 32.05; width of head . 5.87, 3.81, 22,89; length of tailpiece , 
81.70, 59.85, 7.38; length of midpiece, 1.91, 35.62, 37.65. Length of 
tailpiece accounts for most of the phenetic variation. A three- 
dimensional perspective of the projection of the tour OTU’s onto 
the 1 first three components, based on a matrix of correlation 
among the four sperm measurements, is given in Fig. 4, Essen¬ 
tially the same pattern was seen in the distance phenogram. South 
Fork and Guadalupe OTU's were clustered close to one another; 
Beartrap and 4th of July OTU's were loosely grouped. This 
grouping of the populations sampled agreed in part with the 
results from multivariate analyses of cranial morphometric data. 
Small sample size's, however, rendei these results of limited value. 
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Ik;. C— - Tlmr-dimmsionul ptojcuion of the four study population* of Mtrroitts 
jrtt'Xit ami.\ onio tin 1 hist (luce jimik 1 1 >;« I t nmpnnrtiN bused on a rnutnx nl tonche 
(ion unions four spci m mt'aMimm'iiis. Totalities are coded ,»> it) Fij^. I 


Bat ular Morphology 

Ynderson (1960) described the baculum of Microtits mexit anus, 
and I noted no major deviations from his description. Of the 61 
barilla examined. 20 possessed neilhei medial nor lateral pro¬ 
cesses; 52 possessed only the medial process; and 12 possessed 
both. Accordingly, of the five bandar measurements taken, only 
three (bacular length, width ol base, and width of shaft) were con¬ 
sistently available for all spec imens. 

Mean values, followed by range in paremheses and sample si/es 
for selected bandar measurements were (measurements are in 
microns and lo< alities art 1 given in the orde r Beanrap, 4th of July, 
South Fork, and Guadalupe): bacular length, 2.81(2.37-3.76) 18; 
2.98(2.07-3.53) 13; 2.69(2.05-3.21) 15; 2.90(2.52-3.23) 16; width of 
base, 1.35(0.98-1.85) 18; 1.36(0.55-1.87) 13; 1.25(0.71-1.77) 15; 
1.39(0.97-1,71) 16; width of shaft, 0.27(0.16-0.37) 18: 0.28(0.17-0.38) 
13; 0.27(0.17-0.38) 15; 0.30(0.22-0.10) 16. Univariate analysis of 
these three measuremenis revealed no significant differences 
among specimens from the four localities. The large amount of 
individual variation (C Vs, 7.3-30.2) evident in bacular morphol¬ 
ogy, even among individuals of similar si/e (and presumably sim¬ 
ilar age), makes staiistical treatment of bandar measurements dif¬ 
lic nil. Banda enlarge and change somewhat rn shape throughout 
life. Although dividing the bacula into age classes would reduce 
variation due to age, estimating the age of individual microlines 
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is difficult, ami the small number of bacula available rendered an 
attempt to age spec imens unfeasable at this time. 

Although there have been several surveys of microtine bacula, 
tew have treated bacillar measurements statistically. Dearden 
(1958) performed an analysis of variance and reported standard 
errors lor some barular characters in several microtines. His 
results indicated that there were subspec ifir differences in the 
length of the bacular shaft in several subspecies of Lagurus curta - 
tus, although the standard errors were nearly twice those that I 
found in the bacular shaft length of Microtus mexicanus . Coeffi¬ 
cients of variation for bacular characters in the Mexican vole were 
large, compared to the corresponding statistic for any cranial 
measurement. This high degree of variation within a population, 
coupled with small sample sizes, makes interpopulational com¬ 
parisons difficult. As a result, bacular morphology offered no 
information concerning relationships among the' populations 
examined in this study. 

A /orphometric A nalysis 

Fourteen cranial and four external measurements were analyzed 
using the SAS univariate program. Descriptive statistics and 
results of the SS-STP tests for these measurements are given in 
1 able 2. Coefficients of variation obtained for the cranial charac¬ 
ters (3.26-8.95) agu e with (hose values leported by Long (1968) for 
re >dents. 

A two-way analysis of variance was conducted to detect possible 
sexual dimorphism and to determine if differences existed among 
localities for cranial and external characters. Differences were 
found between sexes for depth of skull, total length, and length of 
hind foot below the probability level of 0.05; for interorbital 
breadth, 0.01. Total length of females was consistently greater 
than that of males, and the interorbital breadth was always larger 
in males than it was in females. Within the Beartrap population, 
females had deeper skulls and longer hind feet than did males; in 
the other three populations, males weie larger than females in 
these two characters. 

For diastema length, condylo/ygomatic length, lambdoidal 
breadth, and le ngth of hind foot, there was no significant differ¬ 
ence among localities. Nonoverlapping subsets were found in 
only two characters: the incisive foramen in specimens from the 
Guadalupe Mountains was significantly longer than in those 
from the other three populations; interorbital breadth was larger 


wii m i \i—MicRoi rs mi xicam s 


13 


1'abi.k 2 .—Descriptive statistics derived from two-way analysis of variance for 
external and cranial measurements of Mumms iiK X.it ;mus, (,unip .s of means 
found to be significantly different at VXJM 1 * were tested with the sums of squares 
simultaneous testing procedure to find the nonoverlapping subsets, (iroups of 
means that were found to be not significantly different at P>£105 are marked ns. 
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3 

1 | 

7.3 

5-9-8.0 

.04 

6.19 

■1 

1 

7.3 

0.2-8.1 

.06 

6.24 



Rostral breadth 



1 


1.1 

3,7-13 

.03 

4.80 

2 


4.1 

3.7-15 

.02 

4.44 

1 



4.1 

3.7-16 

,02 

1.78 

3 


4.0 

3,0-18 

.02 

132 
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I ab i k 2. —Continued. 


Intrrnrbital breadth (males 


3.4 

2.9-3.9 

.08 

7.26 

3.4 

3.0-8.6 

.03 

4.76 

3.2 

2,8-3.6 

.01 

6.07 

3.2 

2.7-3.9 

.08 

6.66 

orbital breadth (females) 



3.1 

2.9-3.8 

.04 

7.58 

3.3 

2.9-5.8 

.03 

6.87 

3.2 

2.7-8.7 

.01 

6.56 

3.1 

2.4-3. (5 

.05 

6.76 

Zy gomatn 

breadth 



15.0 

13.3-16.7 

.10 

5.00 

14.7 

13.1-16.0 

.08 

4,57 

14.7 

13.4-10.4 

.06 

1.01 

14.7 

12.8-10.1 

.06 

1.60 

/ 1 re lamb o i d a l b rca d t h 



10.5 

9.5-11.7 

.03 

3.51 

10.1 

9.1-11.2 

.05 

3.58 

10.3 

9.4-11.0 

.04 

3.31 

10.3 

9.5-11.8 

.03 

3.56 

Lam bdoidal breadth 



11.7 

10.6-12.6 

.06 

4.03 

11.7 

10.5-12.6 

.04 

3.90 

11.5 

10.1-12.4 

.04 

3.85 

11.5 

10.2-12.4 

.07 

4.70 

Height of skull 



10.0 

9.2-11.0 

.03 

8.29 

9.9 

9.0-11.2 

.0.5 

4.10 

9.9 

9.2-10.9 

.06 

1.10 

9.8 

9.0-10.6 

.03 

3.46 

epth of brainease (males) 



7.8 

7.3-9.2 

,01 

1.21 

7.7 

7.0-8.8 

.05 

1.75 

7.6 

7.1-8.2 

.07 

4.05 

7.4 

6.3-8.1 

.08 

5.42 

pth of hr aim 

use (femalesj 



7.8 

7.1-8.5 

.04 

3.68 

7.7 

6.9-8.5 

.06 

1.87 

7.6 

6.9-8.6 

.05 

1.64 

7.4 

6.7-8.3 

• 07 

5.13 
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TA b i. t: 2 .—Co n tinned. 


Rostral length 


2 

1 

0.3 

1.9-7.3 

.05 

8.23 

1 

1 

1 I 

6.2 

1.9-7.5 

.07 

8.91 

i 


0. 0 

■1.7-7.3 

.07 

8.95 


1 

5,9 

4.7-7.1 

.01 

7.08 



Total length (molest 



i 


127.1 

111.0-148.0 

2.01 

8.23 

9 


131.2 

107.0-150.0 

1.18 

0.87 

3 

Ml 

131.2 

119.0-H8.0 

1.05 

5.80 

■1 


127,1 

114.0-1 10.0 

1.75 

0.30 



Total length {females} 



1 


130,1 

95.0-151.0 

1.91 

8.90 

2 


131.3 

109.0-173.0 

1.03 

9.17 

3 

IIS 

131.1 

1 10.0 158.0 

1.38 

7.84 

4 


120.2 

111.0 1:52.0 

1,82 

8.09 



/ 

.rngth of tail 



i 

1 

30.0 

22.0-42.0 

.48 

11.34 

2 

1 


20.3 

22.0-58.0 

.12 

15.13 

1 


28.7 

15.0-12.0 

.55 

15.30 

3 


28.3 

18.0-39.0 

.37 

13.58 



Length of hind foot {males. 



J 


17.3 

15.0-20.0 

.21 

6.35 

2 


17.8 

15.0-22.0 

19 

8.15 

:s 

11 s 

17.3 

10.0 19.0 

.1-1 

5.79 

4 


17.7 

15.0-20.0 

.37 

9.53 



Length of hind foot {fern ales} 



1 


17.5 

15.0-20.0 

19 

0.70 

2 


17,2 

12.0-23.0 

.22 

9.85 

3 

ns 

17.2 

1 (>.0-20.0 

.12 

5.30 

1 


17.5 

15.0-21.0 

.29 

9.37 



Length of ear 



4 


12.7 

10.0 18.0 

.25 

14.19 

2 



12.4 

9.0-19.0 

.14 

12.28. 

1 



12 1 

10.0-19.0 

.17 

11.55 

3 1 


12.0 

9.0-15.0 

.10 

8.32 


in individuals from Bfartrap and 4th of July populations than 
those from South Fork and Guadalupe localities. Fht remaining 
characters showed from two to three overlapping subsets but pro¬ 
vided no distinct groupings. 

Characters with high coefficients of variation (length of tail and 
length of ear) and the four characters which displayed sexual 
dimorphism (depth of skull, total length, length of hind foot, and 
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mterorbital hieadth) were eliminated from the subsequent multi- 
variate analyses. The remaining 12 cranial measuiement.s were 
then analyzed using the NT-SYS* multivariate analysis program. 
Correlation and distance matrices were computed and phono¬ 
grams representing the phonetic relationships of the four (VI U s 
were plotted. In addition, a Pearson product-moment correlation 
was computed comparing these two morphometric matrices with 
the corresponding electrophoretic matrices. The correlation coeffi¬ 
cients for the two distance matrices and for the two correlation 
matrices were not significant. The distance phenogram for mor¬ 
phometric characters differed from that produced for sperm mea¬ 
surements presented earlier, with two clusters: Beartrap, 4th of 
July, and South Folk forming the first, and Guadalupe forming 
the second. 

The first three principal components were extracted and plot¬ 
ted, yielding the relationships illustrated in Fig. 5. Phenetic varia¬ 
tion expressed by the first principal component was 54.78 per 
cent; second, 37.75; third, 7.47. The major contributing characters 
for each component and their percentage contribution were: first 
principal component: total length of skull (14.48), palatilar length 
(13.27), condylo/ygomatic length (11.34), zygomatic breadth 
(12.95), and lambdoidnl breadth (12.40); second component: total 
length of skull (27.42) and condylo/ygomatic length (20.62); third 
component: total length of skull (15.55) and prelambdoidal 
breadth (18,74). Thus, characterization of the OTU’s by principal 
components is, to a large degree, dependent on total length of 
skull. The three-dimensional plot arranged the OTU’s into two 
groups that corresponded to those of the distance phenogram. 
Flic Beartrap and 4th of July samples were again grouped much 
closer to each other than to the South Fork and Guadalupe sam¬ 
ples. I he latter two populations are separated somewhat along 
axis II. 

Thus, morphometric analyses yielded no conclusive evidence 
concerning relationships of the four populations studied. lire 
four populations examined generally fell into two groups, corres¬ 
ponding to the same two groups produced by analysis of sperm 
morphology. There is some evidence from the 1 phenogram based 
on morphometric characters and the results of the SS-STP test 
that the Guadalupe population is more 1 distinct than the other 
thiee populations. I lie Guadalupe population has a significantly 
longer incisive foramen. However, for interorbital breadth, the 
Beartrap and 4th of July populations have relatively greater inter- 
orbital breadth than the South Fork and Guadalupe samples. 
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Tttf'xicanu a imio ihr fiisi ihitr principal components based on a matiiv t>t rourl.i- 
linti amont* 12 cranial iiu;tsm rmrnts. L<»( al itics hr- 1 oilc-d as in Fig. 1. 

Piincipal component analyses separated the four populations into 
the same two groups identified in the sperm analysis and in the 
morphometric phenogram. 

Mierotine subspecies generally have been delimited by size dif- 
f cremes as well as qu.ilitiative characteristics* such as pelage 
color. Microtux mexuamis guadalupensis was distinguished by 
Bailey (1902) from M. m* mogollonensis primarih on the basis of 
cranial characters. Bailey (1902) reported measurements for the 
male holoiype of M. m r guadalupensix , six of which may be com¬ 
pared with those taken in this study: toial length, 152; length of 
tail, 31; length of hind foot, 20; basal length of skull, 24.5; length 
of nasals, 7.5; and zygomatic breadth, 16.0. A comparison of these 
moasutemems with those in Table 2 reveals that the total length 
of the holotype exceeds that of any Guadalupe male examined, 
and equals that of the largest Guadalupe female examined. The 
tail of the holotype was longer than the average length of tail 
found in this study, and the length of hind foot equals the largest 
value reported herein I he three skull measurements of the holo¬ 
type all are near the average for those measurements in the Guad¬ 
alupe specimens examined. 


Elect roph oret ic A na Iysis 

Nineteen protein systems were investigated, but only i 6 were 
scored with confidence. These systems contained 24 storable loci, 
which are listed in Table 3 along with the frequencies of each in 
the population. In the total sample, 14 loci (58 per cent) were 
polymorphic, and five (LDH-2, ADH G-G-P, MDH 1. and FS f-7) 
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Cable 3 .—Alleles and frequencies (in parentheses) at 24 loci in Minorus mexita- 
lius and the mean proportion of individuals heterozygous at cadi locus (h ); h is 
ai'eraged freer all four populations (X=S0). 


1 u\ 

Bcailiup 

-lib of Julv 

Snulh I'mk 


h 

LD1I-1 

100 (1.00) 

100 (1.00) 

100 (1.00) 

100 (1.00) 

0.00 

EDI 1-2 

100(.975) 

72(.025) 

100 (1.00) 

100 (1.00) 

100 (1.00) 

0.01 

MDH-1 

100(.975) 
7 (J(.0 23) 

100 (1.00) 

100 (1.00) 

100 (1.00) 

0.01 

MDH-2 

100 (1.00) 

-100(1.00) 

-100(1.00) 

'100(1.00) 

0.00 

ME 

100 (1.00) 

100 (1.00) 

100 (1.00) 

100 (1.00) 

0.00 

IDU 

100 (1.00) 

100 (1.00) 

100(.925) 

L00( .650) 

0.19 




85(.075) 

85(.350) 


SDII 

100 (1.00) 

100 (1.00) 

100(.950) 

116(.050) 

100 (1.00) 

0.03 

ADH 

-100(1.00)) 

—100(.825) 

—100(.950) 

-100(1.00) 

0.11 



—86(.175) 

—86(.050) 



LAP 

100 (1.001 

100 (1.00) 

100 (1.00) 

100 (1.00) 

0.00 

PC a 

IUU{.7f>0) 

100(.800) 

lOO(.HOO) 

100(.250) 

0.40 


54(.250) 

33(.200) 

51(.200) 

54(.750) 


P( ;M 

100(.825) 
160(.175) 

100 (1.00) 

100 (1.00) 

100 (1.00) 

0.06 

VB 

100 (1.00) 

100 (1.00) 

100 (1.00) 

100 (1.00) 

0.00 

IPO 

100 (1.00) 

100 (1.00) 

100 (1.00) 

100 (1.00) 

0.00 

GOT-1 

100 (1.00) 

100 (1.00) 

100 (1.00) 

100 (1.00) 

0.00 

GO 12 

-100(1.00} 

—100(.900) 

— 100( .700) 

“100(1.00) 

0.15 



-41 (TOO) 

—41 (.300) 



G-6-P 

100 (1.00) 

100 (1.00) 

100 (1.00) 

I00(,975) 

8«(.025) 

0.01 

(.DILI 

100 (1.00) 

100 (1.00) 

100 (1.00) 

100 (1.00) 

0.00 

ODIL2 

100 (1.00) 

100 (1.00) 

100 (1.00) 

100 (1.00) 

0.00 

GDH-3 

100 (1.00) 

100 (1.00) 

100 (1.00) 

100 (1.00) 

0.00 

aGPDILl 

100 (1.00.3 

100(.950) 

115(.050) 

100 (1.00) 

100 (1.00) 

0.03 

a(TDII-2 

100(.925) 

10O(.975) 

100(.950) 

100(.975) 

0.09 


129(.050j 

90(.025) 

129( .025) 

129(.050) 

129(.025) 


aGPDH-3 

100(.925) 

100 (1.00) 

100(.950) 

100 (1.00) 

0.04 


N2(.075) 


142(.050) 



ESI -1 

! 00(.-175) 

100(.800) 

100(.850) 

100C.925) 

0.43 


93(.525) 

93 (.200) 

93(. 150) 

93 (.075) 


ESI'-7 

1U0C.975) 

100(.975) 

100 (1.00) 

100 (1.00) 

0.03 


05(.025) 

65(.025) 





were present in a frequency of only five per cent. Acid phospha¬ 
tase (AGP), alkaline phosphatase (AKP), and xanthine dehydroge¬ 
nase (XI)H) were incompletely scored, but did exhibit some 
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polymorphism. The major features of the allo/.yme \ariations 
observed in Microtus tnexicanus are listed below. 

Alalate dehydrogenase (MDH).—Two MDII loti were observed; 
the cathodally migrating MDH-2 was rnononiorphie for all popu¬ 
lations. MDH-1 exhibitt'd two alleles; the MDH-l 70 allele was 
present in a single specimen horn Beartrap. 

Isodtrate dehydrogenase (IDH).—This system was represented 
b\ two alleles. 1 he IDH 85 allele was found in three individuals 
Irom South Folk and in 13 individuals from the Guadalupe sam- 
plr. 

Sorbitol dehydrogenase (SDH).—Two SDH alleles were identi¬ 
fied; the SDH Ilb allele was found in two individuals from South 
Fork. 

lie ohol dehydrogenase (ADH). — Two ADH alleles were 
observed in this cathodally migrating system. The ADH Ho allele 
was present in seven indis iduals from 4th of July and in two 
individuals from South Fork. 

Phosphnglucose isomerase (PGI).—Three alleles were repre¬ 
sented in (his system. The PGI 100 allele otcuned in all samples. 
PGI 4 was identified in eight individuals from Beartrap, eight 
from South Fork, and in all individuals from Guadalupe. PGI 3 * 
was present in eight specimens from the ith of July Campground. 

Phosphoglucomutase (PGM).—This locus was represented by 
two alleles. All populations exhibited the PGM 100 allele; PGM 160 
was identified in six individuals irom Beartrap. 

Glutamate oxalate transaminase (GOT).—Two loci were 
observed; GOT-J was rnononiorphie for all populations, and the 
polymorphic GOT-2 allele migrated cathodally. Four individuals 
from the 4th of July sample and 10 from South Fork exhibited the 
(.01-2 41 allele. 

Glue ose-6-phosphate dehydrogenase (G-b-P).—Two alleles were 
observed; the G-6-P lu " allele was represented in all populations; 
the G-6-P* 8 allele was found in one individual from the Guada¬ 
lupe sample. 

aGlycerophosphatc dehydrogenase (aGPDH). — Three loci were 
identified and all exhibited polymorphism. The aGPDH-l lh 
allele was found in two individuals from ith of July. «GPDH~2 
was rcpiesented by thiee alleles; all populations exhibited the a 
GPDII-2 10 " and aGPDII-2 1 ^ alleles, and the aGPDII-2 90 allele 
was present in only one specimen from Beartrap. aGPDH— 3 142 
was detected in two specimens fiom Beartrap and in two from 
South Fork. 
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I'ahi.k 1 .—Proportion of individuals heterozygous f>er locus per population (hi. P 
is the proportion of 24 Ion polymorphic in each population (loti with t ommunrst 
allele at a frequency of 'XK V5 were considered monomorphic g //, mean proportion 
of loci heterozygous per individual, ,Y = 20. 


1 J K IIS 


Popuhilious 

siimplrt! 



lih ol |ul> 

South J-i >t k 


J .liH 2 

.05 




MD1M 

. 05 




IDH 



.15 

.60 

SDlt 



.10 


\I)H 


.35 

.10 


PGI 

.30 

A0 

10 

.50 

PCM 

.25 




GOl 2 


.20 

10 


G-6-P 




.05 

aGPDH-1 


.10 



cvGPDH-2 

.15 

.05 

.10 

.05 

ftGP DH-3 

.05 


.10 


EST-1 

.-if) 

.10 

.20 

.15 

EST-7 

.05 

.05 



P 

.166 

.208 

.353 

.125 

H 

.051 

.065 

.065 

.056 


Esterases (ESI).—A total of seven esterase systems were identi¬ 
fied, but only two could be scored with confidence. EST-1 con¬ 
tained two alleles found in all four populations. EST-7 also was 
composed of two alleles; EST-7 65 was found in one specimen from 
Bear trap and in one from 1th of July Campground. Of the five 
remaining esterases identified. ES1-2 and ESI-3 exhibited 
polymorphism in all populations. EST-4 and EST-5 were poly¬ 
morphic in all populations except that from the Guadalupe 
Mountains, and EST-6 was monomorphic in all populations. 

Monomorphic proteins. — The following proteins were scored as 
monomorphic in all populations sampled: lactate dehydrogenase 
(L1)H); malate dehydrogenase-2 (MDH-2); malate enzyme (ME); 
leucine arninopeptidase (LAP); albumin (AB); indolphenol oxi¬ 
dase (IPO); glutamate oxalate transaminase-1 (GOT-l); and glu¬ 
tamate dehydrogenase-1, —2, (GDIIT, GDI 1-2, Gl)H-3). 

Genic ^Table 1 gives the proportion of polymor¬ 

phic loci Jii each population, the proportion of heterozygous loci 
per individual, and the proportion of heterozygous individuals 
per locus per population. Genic heterozygosity (H) is based on the 
analysis of 21 loc i (11 polymoprhic. 10 monomorphic). 
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Tabi.k 5.—Coejfifirnts of genie similarity {Rogers' S), getite identity (t\n\s h and 


genetic distance. 

f Xri’s hi 

respectively, for 
nu*\ii .iiiuv 

four populations 

of Mi< itiius 


Hi-nm li[> 

kli >>f |iily 

South Uilk 

Ciumhilujjr 

lit 11 It.l]] 

1.0(100 

0 .9535 

0.0401 

0.0270 


0.0000 

0.0HSM 

0.9847 

0.9520 


0.0000 

0,0117 

0.9155 

o.tMiifi 

■Hh of July 


] .0000 

0.0572 

0.01 12 



0.0000 


0.95-13 



0.0000 

0.0007 

0.0I OH 

South Fotk 



i .000(1 

0.9! 19 




0.0000 

0.0500 




0,04)041 

0.0 1H 

(iwuukiiupe 




1.0000 





0,0000 





0.0000 


The loci that contributed most to tin heterozygosity values var¬ 
ied among populations. Phosphoglucose isomera.se (PCI) was a 
major contributor in all populations, wheieas PST-1 was respon¬ 
sible for most of the individual heterozygosity in the Beat trap and 
hh of July populations. Othet loci connihuting significant 
amounts of heterozygosity were PGM in Real trap, AD11 in 4th of 
July, GOT-2 in South Fork, and II)H in the Guadalupe sample. 

The Bearnap population exhibited lour unique alleles (LDH- 
2 U , MDH-l™, PGM 1 '"’, and aOPDIl-^ 0 ), wheiras tin- -1th nr July 
population possessed only two unique alleles (I*GI and 
cvGPDH 1 1S ). The remaining two populations had one unique 
allele each; South Fork (SOFT 116 ) and Guadalupe (G-6-P 8 *). 

Genic similarity ,—Coefficients of genic simikuity betw een pop¬ 
ulations ware calculated, using Rodgers’ similarity, S (Rodgers, 
1972) and Nei’s identity. / (Nei, 1972). Values for / aie generally 
slightly highet than those tot S\ but both measures give compara¬ 
ble results. Both are reported in Table 5. 

Electrophoretic studies genet ally show high levels of polymor¬ 
phism in natural populations, with reduc ed levels of heterozygos¬ 
ity in small, isolated populations. Selander et at. (1971) demon¬ 
strated low levels of genu vai lability (II 0.018} in insular 
populations of Pernmyscus polionotus compared to the larger 
mainland populations (H = 0.054-0.088). Similar tesults were 
obtained for sev eral species of Perornyscus (Avi.se et at ., 19746), in 
which insular subspecies bad an average of less than one per cent 
heterozygous loc i. Reduced levels of hetc tn/ygosity in small, iso¬ 
lated populations arc 1 thought to be due to genetic drift. 
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The level of genetic variability is relatively uniform for the four 
populations of Mexican voles studied, ranging from 0.051 to 0.065 
(mean, 0.060). These values are consistent with Jreterozygosity 
values for other mainland vertebrate populations (II “ 0.01-0.09) 
as reported by Selander and Johnson (1975). The proportion of 
the 24 loci examined that were polymorphic (P) ranged from 
0.125 to 0.553 (mean, 0.213). This figure also agrees with the value 
of 0.202 that Selander (1976) listed as an average value for rodents. 

Rodgers' coefficient of genic similarity (S) between the Guada¬ 
lupe sample and the other three populations studied is less than 
0.93, whereas those between the other three all are more than 0.91. 
This relationship is more pronounced in Nei’s identity coefficient, 
where the values separating the Guadalupe sample and the 
remaining populations aie less than 0.96; those between the other 
three are greater than 0.98. Coefficients of similarity generally 
range from 0.90 to 1.00 for conspecific populations of todents, 
although in strongly divergent populations of Peromy.scus polio- 
notus on Florida’s barrier islands, the average similarity value 
drops to 0.84 (Selander and Johnson, 1973). 

Other biochemical studies comparing subspecies have yielded 
Rodgers’ similarity values of 0.769 in Mas ?nusculus (Selander ft 
al., 1969), 0.89 to 0.95 in Peromy&cus boylii, and 0.75 to 0.79 in P. 
pohonotus (Avise ft al. , 1974a), and 0.86 in P. eremicus (Avise et 
al., 19746). 1 he low values for AIus in reality could be due to the 
fact that the two subspecies studied are incipient species, whereas 
those reported for P. polionotus result from comparison of a sub¬ 
species occurring on Florida’s barrier islands to a mainland popu¬ 
lation. 'The values obtained in this study ranged from 0.9119 to 
0.9572, and, because these are well wdthin the range of values 
reported for other rodent subspecies, would seem to indicate that 
the populations examined are no more differentiated than would 
be expected on the basis of their current taxonomic rank. Unfor¬ 
tunately, no data from other microtines are available for compari¬ 
son. 

Electrophoretic analysis revealed a close association among all 
lour populations studied, but also indicated that there could be 
some slight differences between the Guadalupe population and 
the remaining three. 

Inconsistencies in the morphometric and electrophoretic anal¬ 
yses result from two different approaches to the same piohiem and 
are to be expected. Electrophoretic analysis frequently fails to dis¬ 
tinguish between or among subspecies that were described 


WILHELM — MICROTI'.S MKXK VXi S 


initiallv on the basis of classical systematic criteria. Whether 
many subspecies are aibiliary units that do not reflect major gene 
differences or whether the resolution resulting from electropho¬ 
retic techniques is insufficient to detect the differences is not 
known at present, hut both elements probably are involved (Avise, 
1975). Contrasts between genic similarity and organismic similar¬ 
ity are not unknown, and it might be a general rule that orga¬ 
nismic evolution and stuctural gene evolution pioceed at vittually 
independent rates, as suggested by Wilson (1976). Therefore*, in 
this study, electrophoretic evidence was accorded less considera¬ 
tion than ware other data in determining systematic placement ot 
the four populations. 


Ch rono logx 

Findley (1969) explained many of the present mammalian dis¬ 
tributional patterns in the southwestern United Slates on the basis 
of a series of boreal expansions and contractions in the late Pleis¬ 
tocene, the comrauions leaving isolated boreal habitats and popu¬ 
lations of certain mammals on scattered mountain ranges. The 
Mexican vole was probablv widespread over the Mexican and 
Colorado plateaus during the cooler, more rnesic pluvial periods 
of the late Pleistocene. The aridity of interpluvial periods caused 
fragmentation of the distribution of that vole, and with increasing 
atidiiy, those fragments became restricted to relatively mesit 
inouiitaintops (see also Sinaitt, 1977). 

The chronolog) of late Pleistocene events in the southwestern 
United States was estimated by Wendorf (1975), based on pollen 
profiles from the Llano Estacado of West Texas. Evidence of a 
continuous boreal forest on the Llano as low' as 1000 meters din¬ 
ing the Eatly Tahoka Pluvial (17,000 BP) corresponds to a depres¬ 
sion of vegetative zones of from 1300 to 1500 meters below their 
current levels. Because inteimontane altitudes in the area under 
study range from about 1250 to 2050 meters, the presence of more 
or less coni in turns boreal forest there in the Early Tahoka seems 
probable. 

The Late Tahoka Pluvial, dated approximately 11,500 BP, evi¬ 
dently was the most recent period of extensive pine-spruce forest 
in this area (Wendorf, 1975). Ibis rather open forest covered 
much of the Southwest, as evidenced by the fact that the Llano 
Estacado was at least 50 per cent covered be spruce-pine forest 
(Wendorf, 1961). During this pluvial period, then, the Mexican 
vole was probably still widely distributed over Arizona, New Mex- 
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ico, and western lex as. M, mcxuanus was present on the Llano 
Estacado at this time in extreme eastern New Mexico, at an eleva¬ 
tion of 1350 meters in what is now Roosevelt County (Slaughter, 
1975). 

Fragmentation in the distribution of A Ucrotus ?ncxiranus prob¬ 
ably began with the Sc harbauer Interval (10,500 BP), which 
marked the end of the Late Ealioka Pluvial. This was a time of 
increasing aridity and was marked by a decline of the pine-spruce 
forest it) the Southwest, although subsequent periods of resur¬ 
gence of coo let and more moist conditions could have allowed 
temporary expansion of boreal habitats. The boieal extensions, 
both in time and geographic scope, since the late Tahokan are 
not presently known, but in all likelihood they were not of suffi¬ 
cient magnitude to connect previously isolated montane habitats. 
The climate of the Southwest has been one of inc reasing aridity 
since the Lubbock Subpluvial (9500 BP) and has undoubtedly 
resulted in the progressive geographic isolation of boreal elements 
on mountain tops. 

Hie subspec ies of Microtus mcxicamis in this study almost cer¬ 
tainly result from 9000 to 10,000 years of isolation (Slaughter 
placed the isolation of M. mexicanus at post-10,000 BP). I he 
degree of differentiation observed between the two nominal taxa 
provides, therefore, some idea of the gross rate of evolution in this' 
species, and, in a general way, underscores the relationship 
be tween geographic isolation and evolution in a small mammal. 

Relative divergence times (T) for paired combinations of all 
populations were calculated following Net’s (1971a) technique 
using the expected number of amino acid differences (D) per pro¬ 
tein that can he detected by electrophoresis. Nei (19715) renamed 
the function D, calling it “genetic distance,” and later presented a 
detailed description (Nei, 1972). Genetic distance is defined as D = 
— logcL where I is Neds identity value. The values of Nei’s genetic 
distance are given in Table* 5. Relative divergence times were then 
estimated using the formula T — Di Di, whe re Di and Eh are the 
genetic distance values for the population pairs in question, and 
are presented in Table 1 6 Using this table, it is possible to postu¬ 
late the order of divergence of the lour populations by listing the 1 
population pairs from one column in order of increasing diver¬ 
gence times. Ibis suggests the following sequence of events. Isola¬ 
tion of the Guadalupe population (4) occurred first. It initially 
was separated from the Beartrap population (1), followed quickly 
h) separation from the 4th of July (2) and South Fork (3) popula- 
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tions. The next population to be isolated was that at Beartrap. It 
was first separated from South Fork, then from 1th of July. The 
last separation isolated the 1th of July and the South Fork popu¬ 
lations. It ran be seen horn the geographic relationships of the 
four study areas (Fig. 1), that the isolation of Guadalupe from 
South Fork probably occurred at about the same time as separa¬ 
tion of the Guadalupe-lth of July and of the Guadalupe-Beartrap 
populations. 

Inasmuch as the calculation of relative time of divergence is 
intimately related to Nei’s identity value (I), it is not surprising 
that these values correlate perfectly with the hypothetical set pie nee 
(derived from observed divergence of characters) of isolation of the 
four populations studied. The highest identity values are found 
between those populations theoretically in contact for the longest 
period of time. Rodgers 1 similarity values (S) indicate the same 
situation, except that they suggest that the Guadalupe population 
was in contact with the Beartrap population for a slightly longer 
period than it was with either the South Fork or the 4th of July 
populations. 

Althougli the calculation of relative divergence times can offer 
only gross estimates of temporal isolation, the sequent e of events 
indicated by these data agrees with what would be expected from 
examination of the geologic evidence. Other recent studies on 
mammals have shown that Nei’s evolutionary divergence time 
does correlate well with both fossil and morphological evidence 
(Nevo et al. y 1974; Zimmerman et ai, 1975; Kilpatrick and Zim¬ 
merman, 1976). As the Pleistocene glaciers retreated northward 
and the climate ol the Southwest bec ame warmer and drier, boreal 
taxa such as the Mexican vole were forced into isolation on mon¬ 
tane refugia. Such restrictions evidently first took place on the 
southernmost Guadalupe Mountains, was followed by isolation of 
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the Bcartrap population west of the Rio Grande on the San Mateo 
Mountains, and ended with isolation of the 4th of July and South 
Fork populations on the Man/ano and Sacramento mountains, 
respectively. It is significant to note that the population thought 
to have been isolated for the longest period of time is the one 
showing the greatest degree of genic and morphological differen¬ 
tiation. 

Results of this study indicate that Minoius mexicanus is an 
evolutionarily conservative microtine. Only slight change evi¬ 
dently has taken place in isolated populations of the species over 
a time span in which several European microlines have developed 
well-marked karyotypic and morphologic differences. No other 
North American microtine has been the subject of a detailed sys¬ 
tematic investigation, however, and the Mexican vole could prove 
to be an exceptionally conservative species. 

TaxoriOrtiic Conelusions 

Subs peci fie differential ion of Minot us mexicanus in the region 
of study undoubtedly has occurred since the isolation of popula¬ 
tions on mountain ranges in the late Pleistocene. 1 he current 
subspecifk Foundries were established by Bailey (1932), with the 
Rio Grande designated as the line of demarcation between M. m , 
magollonensis (occurring in the mountains to the west of the 
river) and M. m. guadalupensis (inhabiting montane areas to the 
east). This arrangement seems logical in a geographic sense, and 
places the Bearnap sample (San Mateo Mountains) in the subspe¬ 
cies mogollonensis and the remaining three samples (4th of July, 
Man/ano Mountains; South Fork, Sacramento Mountains; and 
Guadalupe, Guadalupe Mountains) in the subspecies guadaiupen- 
sis. 

Results of the sperm and morphometric analyses demonstrate, 
however, that the 4th oi July sample bears a closer relationship to 
the Bcartrap sample than it does to either of the other two popu¬ 
lations. Therefore, the 4th of July population is here transferred 
to the subspecies A/, ?n. mogollonensis. The South Fork and 
Guadalupe populations remain refenable to M. m. guadalufrensis. 

There exists, however, a certain degree of divergence between 
the Guadalupe and South Fork samples, as evidenced by both 
c ranial characters and values of electrophoretic similarity. Consid¬ 
ering the isolated nature' o! these two populations and the degree 
of divergence exhibited, I conclude that they constitute an exam¬ 
ple of incipient subspcc iation. The Guadalupe Mountains are 


W’lLMF.LM—MICRO! VS MEXIC'AXLS 


27 


lower m altitude and considerably drier than are the other three 
mountain masses. This results in extremely restricted areas of hab¬ 
itat for Mexican voles and severely limits the population si/e, 1 lie 
Guadalupe sample thus represents a population with several of 
the attributes typically thought to contribute to evolutionary div- 
eigenre. Geographic isolation, small population si/e, occupancy 
of a heterogenous habitat dividing this population into several 
smallei denies, and a harsh or stringent environment relative to 
other conspecific populations, might be reflected in the degree of 
divergence found in the Guadalupe population oi Microtus mexi- 
canus. 
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